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ABSTRACT 

Galaxy structure and morphology is nearly always studied using the light originat- 
ing from stars, however ideally one is interested in measuring structure using the stellar 
mass distribution. Not only does stellar mass trace out the underlying distribution of 
matter it also minimises the effects of star formation and dust on the appearance and 
structure of a galaxy. We present in this paper a study of the stellar mass distributions 
and structures of galaxies at 2; < 1 as found within the GOODS fields. We use pixel 
by pixel K-corrections to construct stellar mass and mass-to-light ratio maps of 560 
galaxies of known morphology at magnitudes z^m < 24. We measure structural and 
size parameters using these stellar mass maps, as well as on ACS BViz band imaging. 
This includes investigating the structural C AS— Gim-M2o parameters and half-light 
radius {Re) for each galaxy. We further identify and examine unusual galaxy types 
with this method, including compact and peculiar ellipticals, and peculiar galaxies in 
some mode of formation. We compare structural parameters and half-light radii in the 
ACS 2:85o-band and stellar mass maps, finding no systematic bias introduced by mea- 
suring galaxy sizes in zsbo- We furthermore investigate relations between structural 
parameters in the ACS BViz bands and stellar mass maps, and compare our result 
to previous morphological studies. Combinations of various parameters in stellar mass 
generally reveal clear separations between early and late type morphologies, but can- 
not easily distinguish between star formation and dynamically disturbed systems. We 
also show that while ellipticals and early-type spirals have fairly constant CAS values 
at z < 1 we find a tendency for late-type spiral and peculiar morphological types to 
have a higher A{M^) at higher redshift. We argue that this, and the large fraction of 
peculiars that appear spiral-like in stellar mass maps, are possible evidence for either 
an active bulge formation in some late-type disks at z < 1 or the presence of minor 
merger events. 
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1 INTRODUCTION 

One of the most intriguing questions in modern astronomy is 
how the universe assembled into the structures we see today. 
A major part of this question is understanding how galaxies 
formed over cosmic time. A popular and rapidly develop- 
ing method of tracing the evolution of galaxies is through 
examining galaxy morphologies and structures over a range 
of redshifts (e.g. Conselice 2003; Conselice, Rajgor & Myers 
2008; Conselice et al. 2009; Lotz et al. 2008; Cassata et al. 
2010; Ricciardelli et al. 2010; Weinzirl et al. 2011). While 
understanding how the morphologies of galaxies evolve, and 
how matter in galaxies is structured, is fundamental, struc- 
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tural studies of galaxies have thus far focused on measuring 
properties in one or more photometric bands, and using this 
as a tracer of evolution. These types of structural analyses 
have always been measured in terms of relative luminosi- 
ties, such that brighter parts of galaxies contribute more to- 
wards their structure and morphology. To understand galax- 
ies more fully however, it is important to study the evolution 
of galaxy structure in terms of their stellar mass distribu- 
tion, as this better reflects the underlying distribution of 
stars in a galaxy. 

Morphological studies have evolved from initial at- 
tempts to describe the range of galaxy forms during the 
early-mid 20th century, towards modern efforts of linking 
the spatial distribution of a galaxy's stars to its formation 
history. In the move away from visual classifications, the 
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goal has been to quantify galaxy structure/morphology in 
a way such that structures can be measured automatically 
and in a reliable way. There are two broad approaches for 
automated classification of galaxies - the parametric and 
non-parametric methods. 

In the parametric approach, the light profiles of galax- 
ies, as seen in the plane of the sky, are compared with pre- 
determined analytic functions. For example, single Sersic 
profiles are fit to an entire galaxy's light profile. Similarly, 
bulge-to-disc (B/D) light ratios are computed by fitting two 
component profiles to a galaxy. These methods are however 
unable to parameterise directly any star formation or merg- 
ing activity which produces random or asymmetric struc- 
tures within galaxies. The "non-parametric" structural ap- 
proaches have no such presumed light profiles implicit in 
their analyses (e.g.. Shade et al. 1995; Abraham et al. 1996; 
Conselice 1997, 2003; Lotz et al. 2004). This is a more nat- 
ural approach towards measuring galaxy structures, as the 
majority of galaxies in the distant Universe have irregular 
structures that are not well fit by parameterised forms (e.g. 
Conselice et al. 2005; Consehce, Rajgor k Myers 2008). Per- 
haps the most successful and straightforward of the non- 
parametric systems is the CAS method, which uses a combi- 
nation of concentration (C), asymmetry (A) and dumpiness 
(S) values to separate galaxy types (Conselice et al. 2000; 
Bershady et al. 2000; Conselice 2003). 

Besides overall structure, the measurement of the size 
evolution of galaxies, as measured through half-light radii, 
also has important implications for their formation histo- 
ries. For example, many studies such as Trujillo et al. (2007) 
and Buitrago et al. (2008) have measured galaxy sizes for 
systems at z > 1 and have found a population of compact 
spheroid galaxies with number densities two orders of mag- 
nitude higher than what we find in the local universe. The 
absence of a significant number of these small sized (as mea- 
sured by half-light radii), high mass galaxies in the local 
Universe suggests that this population has evolved and has 
perhaps merged with other galaxies. However, these stud- 
ies are measured based on the light originating from these 
galaxies, and it is not clear whether sizes would change when 
measured using the distribution of stellar mass rather than 
light. 

All of these methods for measuring the resolved struc- 
tures of high redshift galaxies depend on measurements 
made in one or more photometric bands. As such, quan- 
titative structures are influenced by a combination of ef- 
fects, including: regions of enhanced star formation, irregu- 
lar dust distributions, differing ages and metallicities of stel- 
lar populations and minor/major mergers. The morphology 
and structure of a galaxy also depends strongly upon the 
rest-frame wavelength probed (e.g., Windhorst et al. 2002; 
Taylor-Mager et al. 2007). Young stars, such as OB stars can 
dominate the appearances of galaxies at blue wavelengths, 
thereby giving a biased view of the underlying mass within 
the galaxy. Longer wavelength observations improve the sit- 
uation, but at every wavelength, the light emitted is from a 
mixture of stars at a variety of ages and the dust content. 
More fundamental is the structure of the stellar mass distri- 
bution with a galaxy, as it is more of a direct tracer of the 
underlying potential. Although galaxy structure has been 
measured on rest-frame near-infrared and /-band imaging. 



which traces stellar mass to first order, we directly investi- 
gate the stellar mass images within this paper. 

We use the method outlined in Lanyon- Foster, Con- 
sehce & Merrifield (2007; LCM07) to reconstruct stellar 
mass maps of galaxies within the GOODS fields at 2; < 1. 
We then use this to investigate the evolution of the distri- 
bution of galaxy stellar mass during the last half of the uni- 
verse's history. We use these stellar mass maps to directly 
measure CAS parameters, and the sizes of these galaxies in 
stellar mass over this epoch. We test the assumptions that 
the CAS parameters and sizes can be reliably measured in 
optical light by comparing the parameters for these galaxies 
measured in t/ieB435 , Veoe , ^775 , and zgbo bands and in stellar 
mass. We finally investigate how these various wavelengths 
and stellar mass maps can be used to classify galaxies by 
their formation modes, revealing how these systems are as- 
sembling. 

This paper is organised as follows: in § [2] and § [3] we 
describe the data, sample selection and method, including 
explanations of the K-correction code we use and the CAS 
analysis. In § [H we present our results, discussing the stel- 
lar mass maps themselves, the comparison between galaxy 
size and the CAS parameters in Z85o and mass. We then 
explore the relations between the structural parameters in 
B435, Veoe, ^775, ^850 and stellar mass. Finally, in §[5l we dis- 
cuss our conclusions and comment on future applications. 
Throughout we assume a standard cosmology of Hq = 70 
km s"^ Mpc"\ and = 1 - Qa = 0.3. 



2 DATA AND SAMPLE 
2.1 Data 

The primary source of our data consists of HST/ACS imag- 
ing from the GOODS ACS imaging Treasury Progran|j. 
The observations consist of imaging in the B435 (F435W), 
^606 (F606W), 2775 (F775W) and zgso (F850LP) pass-bands, 
covering the Hubble Deep Field-North (HDF-N) area. The 
central wavelengths of these filters, and their full- width at 
half-maximum, are: F435W (4297, 1038 A), F606W (5907, 
2342 A), F775W (7764, 1528 A), F850L (9445, 1229 A). 
The images have been reduced (using the ACS CALACS 
pipeline), calibrated, stacked and mosaiced by the GOODS 
team (Giavalisco et al. 2004). The field is provided in many 
individual image sections, with the HDF-N data divided into 
17 sections, each of 8192 x 8192 pixels. The total area of 
the GOODS survey is roughly 315 arcmin^, and we utilise 
imaging which was drizzled with a pixel scale of 0.03 arcsec 
pixel" ^, giving an effective PSF size of ^ O.l^^in the z— band. 

The sample consists of galaxies selected in the GOODS- 
N field with z ^ 1 and a F814W magnitude of < 24. The 
sample was also restricted to those galaxies with data in all 
four ACS photometric bands such that we can subsequently 
obtain the most accurate K-corrections and stellar masses 
possible. The final sample consists of 560 objects, which were 
individually cut out of the GOODS-N ACS imaging, and 
handled as separate entities. This is described in detail in 
the following section. 
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Photometric redshifts are available for the whole sam- 
ple (Mobasher et al. 2004), and spectroscopic redshifts are 
available for 404 out of 560 galaxies, as found by the Team 
Keck Redshift Survey (TKRS; Wirth et al. , 2004). When 
spectroscopic redshifts are available they are used, else the 
photometric values are substituted. 



3 METHOD 

3.1 Stellar Mass Maps 

To create stellar mass maps of our galaxies, each pixel in 
each galaxy image is treated individually throughout the 
analysis. We convert fluxes from counts per pixel to appar- 
ent magnitudes per square arc-second for each pixel in every 
image in all four photometric bands, as well as calculating 
their associated errors. Pixels with negative fluxes were as- 
signed values four orders of magnitude smaller than the typ- 
ical flux so that even pixels with low signal-to-noise values 
can be mapped. This however means that our resulting val- 
ues for the CAS parameters, especially asymmetry which 
uses the background light to do a correction, are potentially 
different than when measured using optical light. 

Because our galaxies are at a variety of redshifts, we 
have to carry out fitting to each SED for each pixel to cal- 
culate the stellar mass within each pixel. K-corrections and 
stellar masses are calculated for each pixel of each image us- 
ing the "K-Correct" code of Blanton & Roweis (2007). This 
is similar to previous methods outlined by Lanyon-Foster 
et al. (2007), Bothun (1986), Abraham et al. (1999), and 
Welikala et al. (2008, 2011). 

The code, K-Correc10 naturally handles large datasets, 
and through its interpretation of the data in terms of phys- 
ical stellar population models, outputs results in terms of 
stellar mass and star formation histories (SFHs) of each 
galaxy. The code allows the input of GOODS data and out- 
puts a stellar mass for each object. The code can also be 
modified so as to treat each pixel as a separate input object, 
which we do in this work. 

The K-correction {Kqr{z)) between bandpass used 
to observe a galaxy with apparent magnitude {ttir), at red- 
shift and the desired bandpass Q is defined as (e.g., Oke 
& Sandage 1968): 



mR = Mq+ DM{z) + Kqr{z) - 5log{h) 
where 



DM{z) = 25 + 5log 



h-^Mpc 



(1) 



(2) 



is the bolometric distance modulus calculated from the lu- 
minosity distance, c/l, Mq is the absolute magnitude and h 
= Ho/100 kms-^Mpc-^ 

The K-correct software contains model templates in 
electronic form and an implementation of the method to fit 
data to models. The code uses the stellar population synthe- 
sis models of Bruzual & Chariot (2003) and contains train- 
ing sets of data from GALEX, SDSS, 2MASS, DEEP2 and 
GOODS. The code finds the nonnegative linear combination 
of N template star formation histories that best match the 

^ distributed at ht tp://cosmo.nyu.edu/blanton/kcorrect7] 



observations using a minimum comparison. With the en- 
tire set of galaxy observations available, K-correct also fits 
for the N template SFHs using a nonnegative matrix fac- 
torisation algorithm. K-correct naturally handles data un- 
certainties, missing data, and deals with the complications 
of observing galaxy spectra photometrically, using broad- 
band filters of galaxies at varying redshifts. 

We use a set of 485 spectral templates to fit 
to our galaxy pixel SEDs using Bruzual & Chariot 
(2003) models with the Chabrier (2003) IMF and Padova 
(1994) isochrones. All of the six metallicities available 
(mass fractions of elements heavier than He of Z = 
0.0001, 0.0004, 0.004, 0.008, 0.02 and 0.05) are used. 

Some of the known areas of uncertainty in the stellar 
population models are in the UV and IR regions. In the UV 
light from young or intermediate aged stellar populations 
can dominate the flux at ^ 1500 A. In the near-IR, thermally 
pulsating asymptotic giant branch (TP-AGB) stars domi- 
nate the flux in some intermediate age populations (Maras- 
ton, 2005). We discuss this issue in terms of stellar masses 
in Conselice et al. (2007) who show that the effects of TP- 
AGB stars do not kick in until at redshifts higher than the 
scope of this study {z > 2), and thus are not a concern in 
the present work. 

Results are outputted from K-correct for each input 
pixel, giving rest-frame fluxes at NUV, U, B, V, R, I mag- 
nitudes plus their associated errors, mass-to-light ratios in 
all of these bands and finally the stellar mass in that pixel. 
We then reconstruct the image of the galaxy in stellar mass 
based on these calculations across the galaxy. We use a 
signal-to- noise limit of S/N = 3 per pixel in stellar mass 
for all optical bands, for a reliable calculation. We also in- 
vestigate how the results would change if pixels are summed 
together, finding similar results with the exception of the 
asymmetry which tends to decrease at lower resolutions. 



3.2 Visual Classifications 

We have classified our entire sample visually, using the zgbo 
galaxy images and the object segmentation maps. We have 
derived the classification system by applying the most ap- 
propriate criteria to this particular dataset. The classifica- 
tion scheme is divided into nine categories, ranging from 
compact objects (spherical with little or no apparent en- 
velopes) to obviously merging systems (with evidence of 
tidal streams and other merger attributes). Six objects in 
the catalogue were found to be unresolved and these were 
removed from the sample. We give a description of these 
types and how they were selected below. We often refer to 
these galaxy types throughout this paper, and Appendix B 
gives a description of these types in terms of our measured 
indices. Figures 2-10 show examples of these galaxy types. 



3.2.1 Compact Ellipticals 

The compact ellipticals, denoted as "cE", were identified 
as compact objects that varied smoothly across their radii 
when viewed morphologically. These objects had little or no 
sign of an envelope usually associated with early-type galax- 
ies. The compact Es are generally found at lower redshifts, 
at z < 0.4 with a few between 0.6 < z < 0.8. The dis- 
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18246 



23041 



25035 



32482 



Figure 1. The number of galaxies in each of our sample classifi- 
cations as described in §3.2. Galaxies were visually classified by 
inspection of both zgbO images and segmentation maps. 



tinguishing feature of these compact Es is their small sizes, 
with half-light radii of < 1 kpc. 



36503 



3.2.2 Ellipticals 

The elliptical galaxies are amongst the most popular and im- 
portant galaxy type in our study, as they are typically the 
most massive galaxies in the local universe. One reason for 
this is that elliptical galaxies, and massive galaxies in gen- 
eral, are the test-bed for understanding theories of galaxy 
formation which often have strong predictions for how the 
most massive galaxies should form (e.g., Bertone & Con- 
selice 2009). The ellipticals we see at high redshift are likely 
the progenitors of the massive ellipticals found today. 

Galaxies with the familiar features of nearby elliptical 
galaxies were classed as such, and given the label "E" . 28 E 
galaxies were identified, which at first glance have the ap- 
pearance of early types, but on closer inspection and by vary- 
ing the image contrast had some peculiar features, such as 
multiple nuclei or minor disturbances in the internal struc- 
ture. These objects were classified as peculiar ellipticals, and 
are labelled as pE (Conselice et al. 2007). 



3.2.3 Spirals/Disks 

Objects with a disc-plus-bulge structure were classified as 
early-type spirals ("eS") if the bulge appears larger than the 
disc component, and late-type spirals ("IS") if the disc is 
larger than the central bulge. Edge-on disc galaxies are de- 
noted by "EO" . 

The spiral galaxies are amongst the most interesting for 
this study given that they often contain two major stellar 
population types segregated spatially. Traditionally this is 



Figure 2. Visual displays of the light (in zgbo, left), stellar mass 
(middle) and the distribution of M/L ratios (right) for the com- 
pact elliptical (cE) galaxies. The M/L map is scaled such that 
white corresponds to higher M/L values, and black to lower M/L 
values. This convention is used in the following Figures [3] to \W\ 



seen as an older stellar population making up the bulge, and 
the spiral arms consisting of younger stellar populations. 



3.2.4 Peculiars and Mergers 

Objects that did not fit easily into any of the previously de- 
fined categories, but whose segmentation map showed them 
to be singular, or unconnected to any apparently nearby ob- 
ject in the image, were classified as peculiar, "Pec". These 
are galaxies that could possibly have merged in the recent 
past. 

Galaxies whose morphology also ruled them out from 
any previously defined class but whose segmentation map 
showed them to be connected to, or associated with, at 
least one other object on the image, but without obvi- 
ous merger signatures, were classified as possible mergers 
("pM"). Galaxies with the pM requirements and which 
showed obvious signs of merging, such as tidal tails, were 
classified as mergers ("M"). 

This information is summarised for quick reference in 
Table [U The sample was independently classified four times. 
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Figure 3. Comparing the light (in 2:850 ), stellar mass maps, and 
the distribution of M/L ratios for examples of peculiar elliptical 
(pE) galaxies. The scale and convention is the same as that used 
in Fig. E] 

The mode of the results were taken to be the true classifica- 
tion and where an indecisive split was produced the lowest 
numerically valued type was chosen. Figure [1] shows a his- 
togram of the sample classifications. 

We note that there are a significant number of peculiar 
galaxies in our sample (e.g., Figure 1), much higher than 
what has been found at lower redshifts in previous work in- 
vestigating galaxy morphology (e.g., Consehce et al. 2005). 
These are systems that at lower resolution would be classi- 
fied as disk- like galaxies in many cases. They are not nec- 
essarily merging systems, but are likely normal galaxies in 
some type of formation, the modes of which are one of the 
features we investigate later in this paper. 

We show examples of each of the above types in Figures 
2-10, with five examples shown for each. These figures show 
the image of the galaxy in the 2:85o-band, the stellar mass 
maps, and the mass to light ratio for each galaxy. 

3.3 CAS Analysis 

We use the concentration, asymmetry, dumpiness (CAS) 
parameters to quantitatively measure the structures of our 



Type 


Description 


E 


Elliptical 


cE 


Compact E 


pE 


Peculiar E 


eS 


Early- Type Spiral 


IS 


Late- Type Spiral 


EO 


Edge-on disc 


Pec 


Peculiar 


pM 


Possible merger 


M 


Obvious merger 



Table 1. Descriptions of the classification scheme. 



19929 



21586 



24707 



30976 



36419 




Figure 4. Comparing light, stellar mass maps and M/L for the 
elliptical (E) galaxies. The scale and convention is the same as 
that used in Fig. [2] 



sample, in all available bands, BViz, and on the stellar mass 
maps. We also measure the Gini and M20 parameters, form- 
ing an extensive non-parametric method for measuring the 
structures and morphologies of galaxies in resolved CCD im- 
ages (e.g., Conselice et al. 2000a; Bershady et al. 2000; Con- 
selice et al. 2002; Lotz et al. 2004; Consehce 2003; Lotz et 
al. 2008). The premise for using these parameters is to tap 
into the light distributions of galaxies, which reveal their 
past and present formation modes (Conselice 2003). The re- 
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Figure 5. Comparing light, stellar mass maps and M/L for the 
early- type spiral galaxies (eS). The scale and convention is the 
same as that used in Fig. [2] 



gions into which the traditional Hubble types fall in CAS 
parameter space is well understand from local galaxy com- 
parisons. For example, selecting objects with A > 0.35 finds 
systems that are highly disturbed, and nearly all are major 
galaxy mergers (e.g., Consehce et al. 2000b; Conselice 2003; 
Hernandez- Toledo et al. 2005; Conselice 2006a). A more de- 
tailed analysis of this problem is provided in Appendix A 
for optical light. 

We measure the structural parameters for the GOODS 
sample using the method of Conselice et al. (2008), with 
slight adjustments made for the stellar mass maps, to en- 
able the code to handle the large values of stellar masses per 
pixel, rather than flux. The radius of each individual galaxy 
within the postage stamp image/mass map is measured on 
the stellar mass map, and we define all our indices within the 
Petrosian radii (e.g., Petrosian 1976; Bershady et al. 2000; 
Conselice 2003) . the Petrosian radius has been found to be a 
better are more reliable radius (and reproducible) than the 
isophotal radius (Petrosian 1976). The limits and relation- 
ship to other radii is described in detail in terms of total 
light in a galaxy by Graham et al. (2005). 

Circular apertures are used for measuring our Petrosian 




Figure 6. Comparing light, stellar mass maps and M/L for the 
late- type spiral galaxies (IS). The scale and convention is the same 
as that used in Fig. [2] 

radii and quantitative parameter estimation. The Petrosian 
radius used to measure our parameters is defined by, 

i?Petr = 1.5xr(7y = 0.2), 

where r(?7 = 0.2) is the radius where the surface brightness 
(or stellar mass per unit area) is 20 percent of the surface 
brightness (or stellar mass per unit area) within that radius 
(Bershady et al. 2000). Note that this is a distance indepen- 
dent measurement, given that surface brightness dimming 
effects both measurements of surface brightness in the same 
way. 

Accounting for background light and noise is extremely 
important when measuring structural parameters, especially 
for faint galaxies, and this must also be dealt with for the 
stellar mass maps. The measured parameters in the B435, 
Veoe, '^775 and 2:850 bands and stellar mass images are cor- 
rected as described in Conselice et al. (2008), by considering 
a background area close to the object, and the segmentation 
map of the object. Using a background close to the galaxy 
itself, any problems introduced by objects imaged on a large 
mosaic, with a non-uniform weight map, and the object itself 
being faint compared to the background, are alleviated. We 
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Figure 11. The change in the asymmetries and concentrations measured on our stellar mass maps as a function of redshift due simply to 
distance effects. Show at the top is the change in the asymmetry parameters for a mixture of nearby galaxies of various types: ellipticals, 
spirals and irregulars, while the bottom two panels are for the concentration index. The left hand side for both shows the change for 
the entire sample we simulation with the red dots and errorbars showing the average change and 1 a variation of that change. The right 
hand side shows these simulations divided up between ellipticals (black solid line), spiral galaxies (blue dotted line) and irregulars (red 
dashed line.) 



review below how the CAS parameters are measured. These 
are described in more detail in Bershady et al. (2000), Con- 
selice et al. (2000), Conselice (2003) and Lotz et al. (2008). 

3.3.1 A symmetry 

We measure the asymmetry of a galaxy by taking the orig- 
inal galaxy image (or stellar mass map), rotating it by 180 
degrees about its centre, and then subtracting the two im- 
ages (Conselice 1997), with corrections for background and 
radius (see Conselice et al. 2000a for details). The centre of 
rotation is found using an iterative process, which locates 
the minimum asymmetry. A correction is made within the 
stellar mass maps for this parameter, as it was found that 
uncertainties in the mass-to-light ratio adversely affected the 
asymmetry measurement. We found that approximately 15 
percent of the stellar mass in the asymmetry calculation 
is left over from random fluctuations in M/L. This was 



compensated for including an addition correction for this 
uncertainty through a slight additional subtraction to the 
asymmetry signal. This was done by scaling the asymmetric 
residuals of the rotated and subtracted image by an increase 
of 15%, that is only 85% of the residual is considered. 

We measure the background light in the same way as for 
the light measures of asymmetry, although because we are 
dealing with a conversion to stellar mass, which is not triv- 
ially done for the background, we have some background val- 
ues which are lower than zero. To deal with this, we placed 
all the negative pixels in the stellar mass map background 
to zero. 

The equation for calculating asymmetry is: 

Where Iq denotes the original image pixels, /iso is the im- 
age after rotating by 180° . The background subtraction is 
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Figure 7. Comparing light, stellar mass maps and M/L for the 
edge-on disk galaxy (EO) galaxies. The scale and convention is 
the same as that used in Fig. [2] 



Figure 8. Comparing light, stellar mass maps and M/L for the 
peculiar (Pec) galaxies. The scale and convention is the same as 
that used in Fig. [2] 



made using light (or mass) from a blank sky area, Bq, and 
is minimised using the same process as for the object itself. 
The higher the values of A, the higher the degree of asym- 
metry the galaxy possesses, the most extreme cases usually 
corresponding to merger candidates (Conselice 2003). 



3.3.2 Concentration 

The concentration parameter measures the the intensity of 
light (or stellar mass) contained within a pre-defined central 
region, compared to a larger region towards the edge of the 
visible galaxy. Concentration is most often defined as the 
ratio of the flux contained within circular radii possessing 
20 percent and 80 percent (r2o, tso) of the total galaxy flux. 



C = 5 X log 



\r2o 



(4) 



A higher value of C corresponds to an object where more 
light is contained within the central region. This measure- 
ment has been shown to correlate with the halo, and total 
stellar masses, of nearby galaxies (e.g. Bershady et al. 2000; 
Conselice 2003). 



3. 3. 3 dumpiness 

dumpiness (S) is related to asymmetry, in that it is used 
to measure the amount of light (or mass) in a galaxy that 
exists in discrete, clumpy distributions. In this definition, a 
smooth galaxy contains light at low spatial frequencies, such 
as elliptical galaxies, whereas clumpy systems have most of 
their light contained in high spatial frequencies. Star forming 
galaxies tend to be clumpy in structure, with high S values. 
We measure dumpiness via: 



10 X 



S(/a3,y ■^x,y) 



(5) 



where, the original image Ix,y is blurred to produce a sec- 
ondary image, I^^y. The secondary image is subtracted from 
the original image to create a residual map, showing only 
the high frequency structures contained within the galaxy 
(Conselice 2003). The residuals are quantified by normalis- 
ing, using the total light in the original galaxy image, and 
then subtracting the normalised residual sky. The smoothing 
kernel, a, used is determined from the radius of the galaxy, 
and has the value a = 0.2 • 1.5 x r{r] = 0.2) (Conselice 2003). 
The centres of the galaxies (roughly the inner 1/lOth the 
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Figure 9. Comparing light, stellar mass maps and M/L for the 
possible-merging (pM) galaxies. The scale and convention is the 
same as that used in Fig. [2] 




Figure 10. Comparing light, stellar mass maps and M/L for the 
merging (M) galaxies. The scale and convention is the same as 
that used in Fig. [2] 



Petrosian radius) are removed during this procedure. This 
parameter is furthermore the most difficult to measure and 
is best measured on galaxies imaged with a high S/N ratio, 
and are only of limited use in this paper. Also, the dumpi- 
ness is highly sensitive to resolution or seeing, although the 
ACS imaging we use at these redshifts is within the range 
where values can be properly compared to nearby calibra- 
tion data sets (Conselice 2003). 



3.34 Gini 

The Gini coefficient (G) is defined through the Lorentz curve 
of the light distribution and does not depend on a predefined 
central position, thus distinguishing it from the concentra- 
tion parameter. If G is zero, the galaxy has a uniform surface 
brightness, whereas if a single pixel contains all the flux, G 
is unity. To calculate G efficiently, the pixel flux (or stellar 
mass) values, /i, are sorted into increasing order and the 
following formula is used: 

G= / ^ (2z-n-l)|/.| (6) 
\f\n{n - 1) ^ 



where n is the total number of pixels in the galaxy (Lotz 
et aL2004, 2008). 



3.3.5 M20 

M20 is anti-correlated with concentration such that galax- 
ies with high concentration have low M20 values. M20 is a 
direct tracer of the brightest regions in a galaxy, but does 
not require a pre-determined central coordinate, rather this 
is calculated to be the location of minimal M20 • M20 is more 
sensitive to merger signatures than concentration and is nor- 
malised by the total moment of the galaxy. It is defined as 

M20 = loglO f , while y < 0.2/,ot (7) 

\ Mtot I ^ — ^ 

n n 
i i 

where Xc, yc is the galaxy centre (Lotz et al. 2004). 
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Figure 12. The colour-magnitude relation for our sample of 
galaxies. Displayed is the (Veoe — ^850 ) colour vs. the magnitude 
of the galaxy as observed in the zg^o band. Different galaxy types 
are displayed, including: compact ellipticals (cE), peculiar ellipti- 
cals (pE), ellipticals (E), early- type spirals (eS), late-type spirals 
(IS), edge-on disk galaxies (EO), peculiar galaxies (Pec), possible- 
merger galaxies (pM) and merging systems (M). This key is used 
throughout the remainder of this paper. 
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Figure 13. Colour in observed (Veoe — ^850 ) vs. Redshift (z) for 
our sample. The points displayed are the same as in Fig. [12] 



3.4 Simulations of CAS parameters 

One issue that we must address within this paper is how well 
we are measuring structure on the stellar mass maps for our 
sample of galaxies, and importantly how the measurement of 
stellar mass structure would change with redshift due only 
to cosmological effects. This is critical if we are to make any 
inferences of the evolution of galaxies in terms of their stellar 
mass distribution over time. This issue has been addressed 
before for visual images of galaxies in Conselice (2003) and 
in Beauvais & Bothun (1999) for velocity fields in spiral 
galaxies. Here we address it for the stellar mass maps. 

We simulate how our stellar mass map creation process, 
and CAS measurements would vary as a function of redshift 
due to a decrease in the signal to noise and resolution. We 
carry out this process in the same way we have for the actual 
galaxies in the GOODS fields that we analyse in this paper. 
We take a sample of 82 nearby galaxies of various types - 
ellipticals, spirals, and irregulars, convert their flux values 
to stellar mass. We also go through the entire process we do 
on GOODS, including setting equal to zero the background 
values with negative stellar masses. 

This stellar mass conversion is done after the galaxy is 
redshifted, so that we are mimicking as much as possible 
how the real observations and measurements are done. The 
simulation itself is done by simulating each image at some 
redshift zi to how it would appear at Z2, with Z2 > zi. In our 
case zi 0. When carrying out these simulations of placing 
lower redshift galaxies to high redshifts we calculate first 
the rebinning factor, b, which is the reduction in apparent 
size of a galaxy's image when viewed at higher redshift. The 
other major factor is the relative amounts of flux from the 
sky and galaxy, and the noise produced from the galaxy, sky, 
dark current, and imaging instrument (e.g., read noise) from 
ACS. The process and details for how these simulations are 
done can be found in Conselice (2003). 

In summary, the surface brightness of the simulated im- 
age must be reduced such that the equation. 



Aiiaz^Nz^Pz^(l + zi) 



(9) 



holds, whereby the galaxy as observed in one filter, and is 
simulated in another with central rest-frame wavelengths of 
Xzi and Xz2 and widths AA^^ and AA^-^. In the above equa- 
tion Nz is the total number of pixels within the galaxy at z 
and pz is the average ADU counts per pixel. The calibration 
constant az is in units of erg s~^ cm~^ ADU~^. 

A sky background, and noise from this backgrounds is 
added to these images by Bz x tz, where Bz is background 
flux in units of in units of ADU s~^. For ACS simulations 
we take these values from the measured background based 
on GOODS imaging (Giavalisco et al. 2004) checked to be 
consistent with the values from the ACS handbook. Other 
noise effects are then added, including read- noise scaled for 
the number of read-outs, dark current and photon noise from 
the background. Our resulting images are then smoothed by 
the ACS PSF as generated by Tiny-Tim (Kriss et al. 2001), 
although using PSFs measured from stars give the exact 
same results. 

From these simulated images at redshifts z = 0.5,1 and 
2, we then reconstruct the stellar mass images in the same 
way as we do for our original galaxies, and then measure the 
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CAS parameters on these stellar mass images. The results of 
this are shown in Figure [TTI for the concentration and asym- 
metry parameters. The M20 and Gini parameters however 
have similar behaviours, as do the galaxy half-light radii. 
We find that the dumpiness indices have a similar average 
difference, but with a much larger scatter. 

What we find overall is that the trend with redshift is 
such that the concentration and asymmetry parameters on 
average do not change significantly when measured in stellar 
mass (Figure llip . We show in the right panel of Figure [TT] 
the change in C and A when divided into early/late/peculiar 
types. Again, on average there is not a significant change 
with redshift, although individual galaxies clearly can have 
significant differences at higher redshift. 

We also investigate how our parameters change at each 
redshift between the stellar mass and the visual images after 
the simulation. We find very little difference except that 
at the highest redshifts of our simulations, at z = 2, we 
find that the average change in the asymmetry parameter 
changes 6 A = {Aopt — ^mass) —0.1, such that the stellar 
mass image is more asymmetric than the optical light one. 
We also investigate the same quantities that we use in later 
figures, finding that twice the difference in the asymmetries 
of the optical and stellar mass image, divided by the sum 
of these, increases steadily until it reaches a values of —1.1 
at 2; = 2 similar to what we see later in this paper for the 
actual values. However at 2; = 1 we find that this change is 
less and closer to —0.1. 



4 ANALYSIS 

The following sections describes the analysis of our sam- 
ple, and what we can learn from examining their structures 
in stellar mass maps, and how this evolves over time. We 
first examine the distribution of visual classifications for our 
sample and their global colours. We later describe the stellar 
mass maps we construct for these galaxies, and then finally 
present a structural analysis of these systems based on the 
distribution of their stellar mass. 

4.1 Global Properties of the Sample 

In Figure [12] we show the observed colour magnitude dia- 
gram for our sample. Each morphological class is plotted as 
a different symbol; solid red circles represent the cE's, open 
red circles show the pE galaxies and the E galaxies are repre- 
sented by red stars. Solid blue squares show the eS and open 
blue squares the IS systems. Edge-on disk galaxies are rep- 
resented by black diagonal crosses and solid green triangles 
show the Pec systems. pM and M galaxies are represented 
by open green triangles and open green stars respectively. 
This key is used throughout the rest of the paper in the 
various figures. 

Based on this, we find that our classified early- type 
galaxies are on average redder and brighter than the pecu- 
liars and pM/M systems, while the spirals tend to be bluer 
across a range of brightnesses, but the differences are not so 
clear as would be expected for a sample at \ow-z only. This 
colour-magnitude diagram does not show such an obvious 
morphological sequence as is seen for purely Hubble type 
galaxies in the nearby universe (e.g., Conselice 2006b). 



The relation between (Veoe — Z85o) colour and redshift 
is shown in Figure 1131 The symbols are the same as those 
in Fig. 1121 As can be seen, there is a general trend for all 
types to become redder with increasing redshift, which is 
an expected feature of redshift. The early types are reddest 
across the whole redshift range, followed by the eS and IS 
galaxies, with the Pec/pM/M systems the bluest across all 
redshifts. The compact ellipticals are blue and faint. Note 
however that some pM and eS galaxies are quite red. What 
this reveals is that there is a decoupled relation between 
the colour and morphologies of galaxies not seen in nearby 
galaxies where this correlation is strong. A specific morpho- 
logical type as measured visually at high-z cannot be used 
to predict the colour of the galaxy. 

4.2 Stellar Mass and Mass to Light Ratio Maps 

By calculating stellar mass to light ratios and the stellar 
mass within each pixel of each galaxy image in our sample, 
we reconstruct the image of each galaxy in terms of stellar 
mass and mass-to-light ratio. We do this for each galaxy, us- 
ing the method described in LCM07, and present five repre- 
sentative examples for each classification in Figures [2] to [TOl 
This figures present galaxy images in Z85o (left), stellar mass 
(centre) and {M/L)b (right), which we discuss below. The 
images are grey-scaled such that the darkest pixels represent 
those brightest in zsso, and those that are most massive. In 
the mass-to-light image, the whitest pixels are those with 
the highest (M/L), or redder in colour. We explain below 
how these maps appear for our various types. 

4-2.1 Compact Ellipticals 

Figure [2] shows five examples of the compact elliptical (cE) 
population in the sample. These galaxies appear similar in 
their stellar mass maps as they do in their zgbo images at first 
glance. The mass-to-light maps are not as uniform, and there 
are important quantitative differences between the stellar 
mass and light images. For example, the galaxy 18246 ap- 
pears to have a relatively higher (M/L) in its core compared 
to its outer parts, whilst galaxy 25035 appears to have a low 
(M/L) in its centre, suggesting that it has a bluer core, likely 
due to star formation. Galaxy 36503 contains a high (M/L) 
ring surrounding a blue core. 

4.2.2 Peculiar Ellipticals 

Five examples of the peculiar elliptical (pE) galaxies are 
displayed in Fig. (3] These galaxies appear more diffuse in 
stellar mass than in zsso, in part due to the loss in contrast 
in the stellar mass maps, but also due to their blue inner 
colours (Fig. [3|). As can be seen in the (M/L) ratio maps of 
these galaxies, they have blue cores, and the effect is most 
pronounced in 22932 and 27429. These galaxies also have 
blue overall colours compared to pure E galaxies (Fig. [T3|) . 
These peculiar ellipticals have been seen and studied before 
in papers such as Conselice et al. (2007). They generally 
have a high asymmetry, and are found amongst the most 
massive galaxies in the universe at z ~ 1. As can be seen in 
the M/L maps for these systems (Fig. [3]), these galaxies also 
have a diversity in how young and old stellar populations. 
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including dust, are distributed in these systems compared 
to normal elliptical galaxies. 



4.2.3 Ellipticals 

Figure |4] presents five examples of early- types from our sam- 
ple. These galaxies appear more diffuse and more distributed 
spatially in their stellar mass than in the zsso band, espe- 
cially at large radii. These galaxies have larger (M/L) ra- 
tios in their centres, except for 36419, which is blue, and 
has a more complicated structure in (M/L) than the others. 
Overall, the E galaxies morphologically have very similar 
structures in stellar mass and 2:850. 



4.2.4 Spirals 

Figures [5]and [6]show the zsso, stellar mass and (M/L) maps 
for the early-type, and late-type spirals. As found in LCM07, 
structures within discs, including prominent spiral arms, are 
often (but not always) smoothed out in stellar mass. This 
is true for early types as well as late- type spirals, as can be 
seen for example in galaxy 25465 (Fig. [5]). Inspection of the 
(M/L) maps reveal that the discs of many spirals are bluer 
than their bulges, as expected. 

The mass-to-light ratio maps of the edge on galaxies 
(Fig. (T)) are mostly homogeneous in structure. However, 
the maps of 19280, 39312 and 49722 all show some patches 
of low M/L. Indeed, 19280 appears to be blue across the 
whole image, which would imply that the dust content of 
this galaxy is low, which is unusual for edge-on galaxies. 



4.2.5 Peculiars & Mergers 

Figure [8] shows the F814W images, stellar mass maps and 
M/L ratio maps for the Peculiar galaxies within our sample. 
The Pecuhar galaxies 22690 and 25584 (Fig. [H) appear to 
contain objects near the primary galaxy, but in both cases it 
is the larger galaxies on the right that is the target. Although 
22690 looks disturbed in zsso , the effect of smoothing in the 
stellar mass map is also seen here. This indicates that dis- 
turbed regions of the galaxy are due to star forming regions. 
The Z850 image of 22690 does not have the regular morphol- 
ogy of a nearby spiral galaxy, although the stellar mass map 
has a structure one would expect of such a galaxy, having a 
central bulge surrounded by a smooth disc. The (M/L) map 
also shows a red central region surrounded by bluer pixels. 

The mass-to-light maps of the Peculiars generally show 
that they are blue and also difficult to see in stellar mass, 
due to their low mass to light ratio. This suggests that star 
formation itself might be difficult to trace in stellar mass. 
However, the stellar mass in the pM galaxy 18917 (Fig. [9|) 
traces the light in the galaxy, despite having a low (M/L). 
Note also that the peculiars are often peculiar spirals, in 
the sense that they look like nearly normal spirals, but with 
some peculiar features. Often these peculiars are quite small 
as well, and many of them are very likely spirals in some type 
of formation. 

The merging galaxy system 21839 (Fig. llOp is simi- 
lar in stellar mass to its zgso image, but possesses a more 
intricate structure in {M/L). The mass-to- light map shows 
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the smaller merging object to be blue, whereas the brighter 
galaxy appears more red but with blue regions in its core. 

The system 29800 (Fig. llOp . whilst clearly appearing 
to be a merger between two galaxies of approximately equal 
brightness in zsso , would likely be classified as a disc galaxy 
in stellar mass. The bottom of the two bright nuclei has 
disappeared completely in the stellar mass image, although 
it is visible as a low (M/L) patch in the mass-to- light map. 

4.3 Comparison Between Galaxy Properties in 
Mass and Light 

In this section we compare how the distribution of stellar 
mass in a galaxy compares to the distribution of light as seen 
in the zsso ACS imaging. One of our goals is to determine 
how appropriate studies in zsso and similar red bands are 
for measuring the mass content and structure of a galaxy, 
and how stellar mass quantitative morphologies differ from 
those measured in light. 

In this section we investigate the relationship between 
galaxy properties in stellar mass and zsso-band light. For 
galaxy size (as measured by the half-light radius) and the 
CAS parameters we plot the normalised difference between 
^850 and stellar mass versus the mean value of zsso and stel- 
lar mass as a representative figure to determine how these 
various quantities change. 

4.3.1 Galaxy Half-light and Half-mass Radii 

One of the basic features of a galaxy is its size, which in this 
paper we quantitatively mean the half-light radii. Galaxy 
half-light radii are found to strongly evolve with time, such 
that galaxies at higher redshifts have a more compact struc- 
ture (e.g., Ferguson et al. 2004; Trujillo et al. 2007; Buitrago 
et al. 2008; Carrasco et al. 2010). However, every size mea- 
surement has been carried out through measurements of 
light, and it is desirable to determine how the half-mass 
radii of galaxies measured in stellar mass maps compares 
with half-light radii measured using light. 

As such, we have calculated half-light (or half-mass) 
radii (Re), in kpc, for our sample in both stellar mass and 
^850 light, and investigate whether half- light /mass radii are 
comparable. We define the normalised size difference as 

2[i?e(^850) - Re{M,)]/[Re{z850) + i?e(M*)], 

and the average size as 

^[Re{z850) + Re{M,)]. 

We use these two calculations so as to avoid biasing the anal- 
ysis when the values become very small in either the stellar 
mass image, or in the zsso band. We plot the relation be- 
tween these values in Figure [Ml The horizontal dashed line 
marks the position of equal size (i?e(z85o) = i?e(M*)). The 
symbols are the same as in Figure [TJ) and the average disper- 
sion for the sample is represented by the black square points 
on the right of the plot. The larger black square points in the 
figure show the average values and the measured dispersion 
for the ranges of | {Re{z85o) + Re{M^)) between zero and 
one kpc, one and two kpc, up to six kpc. This convention 
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i[Re(z) + Re(m)]/Kpc 

Figure 14. Normalised effective-radius difference vs. mean ef- 
fective radius of the sample in zgbo and stellar mass. The error 
bar to the right of the plot shows the average dispersion for the 
whole sample. The large black squares denote the average values 
in equally spaced bins of mean size, with error bars showing the 
dispersions of these values. This convention is used in all of the 
parameter comparison plots. The symbols plotted here are the 
same as in Figure 12. 



is used for each of the parameter comparisons in Sections 

We note that there is a steady progression for increas- 
ing average galaxy size with morphology, from the compact 
ellipticals, early types and early spirals to the late- type spi- 
rals, with the peculiars and edge-on galaxies being more ran- 
domly distributed. The points are scattered about the line 
of equality such that nearly half (42 percent) of the galaxies 
have i^e(^85o) > Re{M^). There is, therefore, a slight ten- 
dency for sizes in masses to be higher, but this does not vary 
significantly with average galaxy size and type (see Table 2). 

There is no clear tendency for any particular morpho- 
logical type to be larger in either the zq^q or stellar mass 
image, which suggests that there is no bias introduced by 
using measurements in zsso band data for size measurements 
(Trujillo et al. 2007; Buitrago et al. 2008). Table shows 
that the average difference between the sizes in the stel- 
lar mass maps and the ^sso-band image is essentially zero, 
demonstrating that the measurements of half-radii in light 
does not differ significantly from the measurements in the 
distribution of stellar mass. 

We also examine the normalised difference in size 
against (Veoe - ^^sso) colour to test whether there is a ten- 
dency for bluer galaxies to have larger radii in stellar mass 
than in light. This is what we might expect to find if star 
formation dominated the light near the centre of the galaxy. 
While, on average, colour does not change with size differ- 
ence, we find that extreme half-radii differences between stel- 
lar mass and zsso are mostly in blue systems. 



Type 


Mean Size 
(kpc) 


Normalised Size Difference 


cE 


0.6 ± (0.3) 


0.14 ± (0.24) 


pE 


1.8 ± (0.2) 


-0.11 ± (0.11) 


E 


1.8 ± (0.3) 


0.05 ± (0.11) 


eS 


2.6 ± (0.4) 


0.00 ± (0.12) 


IS 


3.7 ± (0.6) 


-0.10 ± (0.12) 


EO 


2.7 ± (0.5) 


0.02 ± (0.13) 


Pec 


2.6 ± (0.5) 


-0.09 ± (0.13) 


pM 


2.7 ± (0.6) 


-0.07 ± (0.13) 


M 


3.0 ± (0.6) 


-0.08 ± (0.03) 



Table 2. Comparisons between mean zgso -stellar mass half-mass 
radius {Re) and normalised size difference in zsso and stellar 
mass, organised by morphological type. 



4.3.2 Concentration 

Analogous to the galaxy size comparison in the previous sec- 
tion, we compare the concentration parameter in zsso and 
stellar mass images, as plotted in Figure [15] The dashed line 
shows C(^85o) = C(M*) and the black square points repre- 
sent the average values in equally spaced bins of average 
concentration. 

The average concentration, in both ^sso and stellar 
mass, is generally high for the early types, and lower in the 
late-types and peculiars, with the early-type spirals being 
spread between low and high values. The compact elliptical 
galaxies have a low average concentration, due to the light 
being spread evenly across a galaxy's pixels. 

The average trend shows that the ratio is slightly 
lower in the stellar mass maps for the early types compared 
to the ^850 band concentration. We find that in general Rso is 
lower and R20 is higher, compared to those in the 2:850 -band 
image, although the effect is mostly due to a smaller Rso 
although this effects is quite small. It is possible that stellar 
mass is less centrally concentrated towards the centre of the 
galaxy in the early- types (raising R20), and/or the stellar 
mass is more diffuse at larger radii than the distribution 
of light (lowering Rso)- It is also the case that for many 
galaxies more star formation is seen in the outer regions 
of the galaxies in the zsso-band, but these pixels are less 
dominant in stellar mass, raising the value of Rso creating 
higher values of C as seen. 



4.3.3 Asymmetry 

Before we discuss the comparison of our stellar mass map 
asymmetry values to the ^sso-band asymmetries, we note a 
few things. First, the lowest asymmetry values for the sample 
are negative in all bands, mainly due to the background 
correction which for very symmetric galaxies will sometimes 
be larger than the asymmetry in the galaxy itself. This has 
the effect, as seen in Figure [161 of skewing the low average 
asymmetry values differences for the smallest stellar mass 
measured asymmetries. As can be seen in Figure 16, except 
for these galaxies with low asymmetry values, there is a clear 
tendency for galaxies to be more asymmetric in stellar mass 
than in zsso, this being the case for a large fraction of the 
non-early type sample. The stellar mass maps also have a 
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Figure 15. Comparison between the concentration parameter in 
the Z850 band and within stellar mass maps. The error bar to the 
right of the plot shows the average dispersion for the whole sam- 
ple. The large black squares denote the average values in equally 
spaced bins of mean concentration, plus dispersions. The symbols 
used here are the same as in Figure 12. 

greater spread in asymmetry values than zsso, and this can 
be more clearly seen in § 14.4.11 

However, we also note that the blue regions of galaxies 
are not so well traced by the stellar mass maps, and often 
led to difficulties with the image contrast, making the galaxy 
features harder to see, as in the case of the late-spiral 19535 
(Fig. [6]). Further galaxies, such as 38722 (Figure 6), are 
clearly more asymmetric and lopsided when viewed in the 
stellar mass band. In this case, it appears that one spiral 
arm remains while the other disappears. 

There are several reasons why the asymmetry value in 
the stellar mass maps is higher than in the images. First, 
when calculating A, noise tends to be magnified due to the 
process of subtracting images. The simulations discussed in 
§3.4 show that there is a tendency for the average galaxies to 
become more asymmetric in stellar mass measurements than 
in light. This can explain part of this, but we find that the 
average difference of around —1 in the relative difference is 
too high to be accounted for solely by these redshift effects. 
We discuss some of the reasons for why the asymmetries will 
be higher in the stellar mass maps than in light. 

When calculating stellar masses for each pixel we as- 
sume that the M/L ratio is approximately constant in the 
surrounding pixels. This should be the case for early type 
galaxies, especially, due to their uniform stellar populations. 
To test how much this variance in M/L could be affecting 
the A values, we have measured the M/L pixel variance in a 
typical early type galaxy (30976, Fig. The mean M/L for 
a pixel in this galaxy is {M/L)b — 1-29 with a typical stan- 
dard deviation in the surrounding pixels of (Tm/l — 0.27. 
Although a correction has been applied to minimise this ef- 



Figure 16. Comparison between the asymmetry parameter in 
^850 and stellar mass. The points and conventions are the same 
as in Figure 12. Note that the asymmetry in the stellar mass is 
lower than the z— band for only the early types. This is due to the 
method of measuring the asymmetry, where the sky background 
is handled in a different way than for imaging, resulting in higher 
values for later type galaxies (§4.3.3). 

feet (see §3.3.1), it is not enough to account for the high 
asymmetry signal in the spiral galaxy stellar mass maps, 
and thus this asymmetry is likely a real effect due to the 
nature of the calculation of this parameter, and especially 
the background, as described in §3.3.1. 

To further understand why late-type spiral galaxies are 
more asymmetric in stellar mass than ^sso, we have exam- 
ined these galaxies in detail and present four examples of 
the late- type sample in both zsso (Fig. [T7|) . and stellar mass 
maps. Two of these galaxies have low asymmetries in both 
stellar mass and ^sso (25751 and 21448), and galaxy 21448 
in particular, demonstrates the smoothing of spiral arms, 
which we observed previously in nearby galaxies (LCM08). 
However, two of these galaxies have higher asymmetries in 
stellar mass than in zsso (19535 and 34946). Object 19535, 
especially, shows little relation in stellar mass to the zsso 
image. The star forming knots and central bulge, which can 
be seen clearly in zsso, translate to large scale asymmetries 
in stellar mass. There are several possible explanations for 
this effect. The star forming regions in 19535 could truly be 
more massive per pixel than the surrounding disk, due to 
the high masses and densities of gas and dust required to 
form such massive star forming regions. This would lead to 
such regions not being smoothed out in stellar mass, causing 
higher A(M*) values. It is also possible that what we have 
interpreted as star forming knots in the disk in some cases 
could be minor merger signatures, that is, low-mass galaxies 
falling into the object itself. 

There is also a trend with morphological-type with in- 
creasing asymmetry such that early types have low asym- 
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Figure 18. Comparison between the Gini parameter in zgso and 
stellar mass. The points and conventions are the same as in Fig- 
ure 12. 



metry and late-types higher asymmetry. Our findings con- 
firm previous studies who found similar trends (e.g. Con- 
selice et al. 2000). We have also investigated the rela- 
tion between asymmetry in stellar mass and (Veoe — ^sso) 
colour, and examine this after splitting into two redshift 
bands: z < 0.6 and 0.6 ^ z < 1. There is a clear 
trend for galaxies with a higher degree of asymmetry to 
be blue, with ((Veoe - ^85o)|a(m,)<o.35) = 1.01(±0.20) and 
((V606-285o)|a(m>,)^o.35) = 1 .22(di0.27) . This trend has also 
been found for asymmetries measured in light (Conselice 
et al. 2003). 



4.3.4 dumpiness 

The dumpiness values in stellar mass maps have a higher 
degree of uncertainty, as discussed in Section 3.3.3. Overall, 
we find a trend for dumpiness to be greater in stellar mass 
than Z850, and this is the case for 91 percent of the sample. 
However for larger values of average there is a trend such 
that S{M^) S{z85o). 

The dumpiness parameter is similar to asymmetry, but 
picks out small scale features, such as compact star clus- 
ters, as opposed to large scale asymmetries. It is, therefore, 
not surprising that uncertainties in M/L cause a large num- 
ber of the sample to have 5'(M*) > 5'(85o), as variations in 
M/L from pixel to pixel filter through in the calculation of 
S{M^). Such variations affect A{M^) and S{M^) more than 
the other parameters, as these calculations involve the sub- 
traction of images, which amplifies the M/L variations and 
large variations of level within the background. 




i[M,„(z) + U^Jin)] 

Figure 19. Comparison between the M20 parameter in zgbo and 
stellar mass. The points and conventions are the same as in Fig- 
ure 12. Because of the inverse nature of M20, galaxies with posi- 
tive differences on the y-axis are those with more 'concentrated' 
light profiles than in the stellar mass maps. 
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Figure 21. The asymmetry-concentration plane for the stellar 
mass maps. The plotting convention is the same as for Figure 12. 
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Figure 17. Figure showing the 2;850-band imaging (left) and stellar mass maps (right) for four late-type spiral galaxies. Each galaxy is 
labelled with its ID number and A(z85o) value. Each image is 4.5 arcseconds on each side. Shown at the top of each image is the ID and 
the asymmetry parameter measured on the stellar mass maps to be compared with the same value computed in the zs5o band. 



4-3.5 Gini Index 

Figure [18] shows the comparison between the values of the 
Gini parameter in stellar mass and the zsso band. Although 
the early-types are more strongly clustered around equality, 
G{M^) = G{z85o), there is a clear trend for the sample to 
have higher G in stellar mass, with G{M^) > G{z85o) in non- 
early types within our sample. This comparison of G{z85o) 
and G{M^) displays a greater separation between types. The 
late-type spirals show an approximately even spread in Gini 
values, from 0.45 to 0.93, but all of these have a higher 
difference between zgbo and stellar mass than the early- types 
and compact ellipticals. 

Figure 1181 reveals that the Gini index is higher within 
the stellar mass maps than in the zsso band, except for early- 
type galaxies, indicating that most of the stellar mass in 
later morphological types is contained within fewer pixels 
within the stellar mass image. This is a significant difference 
between the early and late- types in our sample. The reason 
for this is that the bright blue regions of these galaxies van- 
ish as the M/L ratio is inversely proportional to L. This 
creates effectively a mass distribution with a larger fraction 
of the mass contained within fewer pixels. For example the 
bulges become more prominent for the late-type disks while 
the arms vanish. The bulges in these systems become more 
prominent and this will rise the value of the Gini index. 

What can often been seen in stellar mass maps of late- 
types and mergers/peculiars is that these bright outer re- 
gions vanish, leaving only the central part of the galaxy 
including most objects in Figure 6 which loose their outer 
parts and appear almost as early- types. Figure 18 also shows 
that the Gini index for the early-types in stellar mass and 



light are similar, likely because of the small variation in the 
M/L ratios of the various pixels. 



4.3.6 M20 

We show the differences between the M20 index in the stel- 
lar mass maps and the zsbo band in Figure [T9l Due to the 
reversed parity in the M20 values, points which lie above the 
dashed line in Figure [TOlhave M2o(M*) > M2o('2^85o), that is, 
the stellar mass is more more diffuse than the light for most 
systems. This would change however if we used a smaller 
radius than the standard definition. Overall, M2o(M*) is 
higher (less negative) than M2o(^85o) for most of the sample. 

The M20 parameter traces the brightest 20 percent of 
the flux in the galaxy, or the greatest 20 percent of the stellar 
mass. This parameter is heavily weighted by the spatial dis- 
tribution of the most luminous, or most massive, pixels, but 
is normalised to remove any dependence on galaxy size and 
total flux/stellar mass. M20, like C, is also not dependent 
on a fixed, pre-determined centre. 

Therefore, if stellar mass exactly follows light within a 
galaxy we would expect that M2o(M*) = M2o(^85o)- Figure 
[19] shows that, while this is approximately the case for early- 
type galaxies and early spirals, the same does not apply for 
peculiars, late- type spirals and edge-on discs. The compact 
ellipticals especially differ from this trend, having M2o(M*) 
< M2o('^85o)- This is certainly due to the fact that these 
compact ellipticals have blue cores which render their stellar 
mass maps more compact in the centre than for galaxies with 
redder cores. 



Galaxy Structure in Stellar Mass and Light at z < 1.0 17 





Figure 20. Asymmetry vs. concentration in the B435, Veoe, ^775 and zgso bands. Solid lines show the classification criteria of Conselice 
(2003). The point in the bottom left segment of the plot shows the average error bar for the whole sample. 



4.4 Stellar Mass Structure 

Concentration vs. Asymmetry 

We examine how galaxies fall in the classic concentration- 
asymmetry plane (e.g., Conselice et al. 2000) to determine 
if different galaxy types, as determine visually can be better 
separated at different wavelengths in this parameter space. 

We plot the relation between asymmetry and concen- 
tration in light for the sample in Fig. [20] and in stellar mass 
in Fig. 1211 Over-plotted are the classification criteria of 
Conselice (2003) and Bershady et al. (2000). Galaxies with 
A > 0.35 are classified as mergers, while those which lie 
above the top line are early- types. Galaxies to the left of the 
middle line in Figure 20 & 21, are labelled mid- types and 
those to right of this line are classified as late- types. 

These criteria for identifying mergers and late-types ap- 
pear more appropriate for our sample in the 2775 and zsso 
bands, while galaxies are more mixed in these regions in 
stellar mass, having many more galaxies with high asymme- 



tries, and spirals are not easily distinguished from merging 
systems. However, the uniqueness of the early- type region is 
more obvious in stellar mass maps, with most E galaxies in 
the 2775 and zsso band plots falling into the mid-type region 
(Fig. EU). 

We find that for our sample, ^sso is the best band in 
which to measure these parameters for separating morpho- 
logical types from each other. It can also be seen in Fig. 
[20] that the A-C relation breaks down when viewed in bluer 
bands, as also shown for nearby galaxies in Taylor- Mager et 
al. (2007). We find that the scatter in A-C increases for all 
galaxy types towards bluer wavelengths, the effect is more 
pronounced for late- type morphologies. 

Figure 1211 shows the A-C relation for the stellar mass 
maps. There is a greater range in asymmetry values in stel- 
lar mass and the criterion suggested by Conselice (2003) for 
merging systems, of A > 0.35, even includes a few E and 
some eS galaxies. There is, however, still a clear separation 
in types between early and peculiar type systems, but the 
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Figure 22. Gini-M2o for the B435, Veoe, ^775 and zgbo bands. The sohd hnes mark the classification criteria of Lotz et al. (2004; Equations 
[I0]and [m. 



spirals (especially IS galaxies) significantly overlap between 
the two. This makes it difficult to distinguish between mor- 
phological types in stellar mass within A — C. However, given 
the fact that many of these spirals appear to have some kind 
of merger or formation mode, the stellar mass structure is 
superior for finding galaxies in active evolution. 

4.4.2 Gini vs. M20 

We plot the relation between G and M20 for light and stellar 
mass in Figures [22] and 1231 The solid lines in Fig 1221 are 
taken from Lotz et al. (2008), who used these parameters to 
classify their sample of normal galaxies and ULIRGs. Equa- 
tion [To] describes the line above which Lotz et al. classify 
their sample as mergers, and Eq. [TT] describes the line sep- 
arating early and late- type galaxies. 

G^ -0.14 •M20 + 0.33 (10) 
G = 0.14 -^^20 + 0.80 (11) 



We again see that the cleanest, in terms of separat- 
ing morphological types, optical relation is in zsso, al- 
though our sample is not best represented by the Lotz 
et al. classifications in Gini-M2o space, with many early-type 
galaxies falling into the merging region of the plot. This is 
partially due to the fact that G and M20 are calculated dif- 
ferently here then in Lotz et al. , as described in §3.3 and 
Lisker et al. (2008). However, we cannot rule out that the 
classifications for this sample may be intrinsically different, 
and it is worrying that all of the edge-on galaxies fall into 
the merging region of G-M20. We also use a similar proce- 
dure to look at the same rest-frame wavelength as Lotz et 
al. does. 

We modify the Lotz et al. (2008) relations for zsso (top) 
and stellar mass (bottom) and show this in Fig. 1231 The 
solid lines in the ^sso plot are our revisions to the Lotz 
et al. relations, based on our sample (Equations [T2]for merg- 
ers and [13] for normal systems). 
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Figure 23. Gini-M2o for 2^850 (top) and the stellar mass maps 
(bottom). In zgbo (top) the dashed lines are the classification 
criteria of Lotz et al. (2004; Equations llOland [TT]) while the solid 
lines mark our revised criteria in zgbo (Equations fT2] and fT3)) . 
In stellar mass (bottom) the Lotz et al. criteria are marked by 
the dashed lines, our zgbo criteria are shown by the dot-dashed 
lines and the solid line shows our separation of early and late 
types/mergers, in stellar mass. 

^(^850) = -0.14 • M20(^850) + 0.38 (12) 
^(^850) = 0.14 • M20(^850) + 0.74 (13) 

while the dashed lines are the original Lotz et al. relations. In 
the bottom panel of Fig. I23[ the solid line shows our revision 
of the merger defining line for the stellar mass maps, (Eq. 

G{M.) = -0.14 • M2o(M*) + 0.45, (14) 

while here our relations for zsdo and those of Lotz et al. are 
shown as dot-dashed and dashed, respectively. 
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Figure 24. Asymmetry in stellar mass maps vs. redshift. The 
solid line shows the fit to the E galaxies (Eq. I15|) , which is approx- 
imately constant across redshifts (z). The dashed line illustrates 
the fit to the late- type spiral galaxies (Eq. [Tg]) . The late- type 
spiral galaxies show a larger A(M^) with increasingly higher red- 
shifts. 



As with the relation for A(M*)-(:7(M*), the G{M^)- 
M2o(M*) relation shows the same general trend as its Z850 
counterpart, but with a larger spread in values. The spi- 
rals, again, are not as neatly defined as in ^850, with the 
eS galaxies separated from the early-types and the late-type 
spirals and edge-on disk galaxies occupying the merging re- 
gion. However, there is a clear separation in the stellar mass 
Gini vs. M20 which does not exist for the optical light. Sev- 
eral of the pE galaxies lie in the merger region, even as 
defined by the higher stellar mass line (Eq. I14p . This has 
also been noted by Conselice et al. (2007), and is likely due 
to these objects having multiple nuclei. 

The G(M*)-M2o(M*) relation shows a clear early/late- 
type split, defined by Eq. [TH as the Sb/Sc/Irr region, as de- 
scribed by Lotz et al. is not appropriate to apply to our sam- 
ple in stellar mass. We discuss other correlations: asymmetry 
vs. dumpiness, concentration vs. M20, size vs. concentration, 
and asymmetry vs. M20 in the appendix A. We evaluate in 
this appendix other scaling relationship between these non- 
parametric relationships, including asymmetry-dumpiness, 
concentration and M20, size and concentration, and asym- 
metry and M20. 



4.5 The Evolution of Galaxy Stellar Mass 
Structure with Redshift 

In this section we investigate the changes in our galaxy stel- 
lar mass maps with redshifts. We take a quantitative ap- 
proach here and find that at z < 1 there is little change 
in the stellar mass structure with redshift, with the excep- 
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tion of the asymmetries of these galaxies which decreases for 
some galaxy types at z < 1. 

We plot A(M*) against redshift (Figure I24p where we 
can see a general trend for galaxies to become more asym- 
metric in stellar mass at higher redshifts. We find that not 
all galaxy types follow this pattern however, with the most 
notable exception being the ellipticals. As shown in Fig. [24] 
there is no trend for the ellipticals to become more asym- 
metric in stellar mass at higher redshifts. 

We note that while the Es remain at approximately the 
same low asymmetries across the range in redshift, there is a 
trend for spiral galaxies to become increasingly asymmetric 
at higher redshifts. We measure this trend by fitting both the 
ellipticals (solid line) and IS (dashed line) galaxies, plotted 
in Figure [23] and quantified in Equations [15] and [16] below. 

Ae{M.) = 0.049(±0.079) • z + 0.021(±0.056) (15) 

is the best fit for the ellipticals and for the late- type spirals, 
the best fit is, 

Ais{M,) = 0.327(±0.093) • z + 0.175(=b0.064). (16) 

Equation [15] shows that the early- type galaxies retain ap- 
proximately constant asymmetries across < z < 1. Equa- 
tion [T6]indicates that there is a relation for late-type spirals, 
such that A{M^) increases with redshift. 

We have already discussed in §4.3.3 how the late- type 
spiral galaxies have stellar mass asymmetries larger than 
their asymmetries in optical light. Since the late- type spirals 
are galaxies which have arms/disks that are larger/brighter 
than the bulge component, these galaxies are therefore dom- 
inated in terms of their light by their spiral structures and 
arms. These galaxies show the most diversity and contract 
between spatial distributions within their M/L and stellar 
mass maps (e.g.. Fig. [6]). 

We note that not all of the late-type spirals have high 
asymmetries in the stellar mass maps, some are as asym- 
metric or less than in the ^sso-band (Fig. [T6|) . Furthermore, 
we find that in general the disk galaxies with higher stellar 
mass map asymmetries have a bluer colour, which is the case 
at both high and low redshifts. In many of these cases, the 
structures are lopsided and/or contain outer features that 
remain after normalising by the M/L map. This interpre- 
tation of stellar mass being more localised than the light is 
consistent with a higher Gini index in the stellar mass bands 
than in the light. 

Some of these asymmetric spirals, as well as many of 
the peculiars resemble the so-called "clump clusters" and 
clumpy spirals found by Elmegreen et al. (2005, 2007), es- 
pecially in stellar mass maps. Their low concentration in 
mass is also similar to the Luminous Diffuse Objects (LDOs) 
found at similar redshifts (Conselice et al. 2004). Indeed, 
there are some late-type spirals galaxies that we would have 
classified differently, perhaps as Pec or pM systems, had we 
performed the classification in stellar mass, rather than in 
Z850- The clumps in galaxies found by Elmegreen et al. are 
also massive, with typical values in the region ~ 10^ — lO^M© 
(Elmegreen & Elmegreen 2005), and reside in young disks at 
high redshift. It is likely that some of our highly asymmetric 
late- type spiral galaxies are similar, or in the same class as 
these clumpy spirals. These features may also be the result 
of minor merger events. 



5 CONCLUSIONS 

We conduct a pixel by pixel study of the structures of 560 
galaxies found in the GOODS-N field ai z < 1 within stellar 
mass maps and Hubble Space Telescope ACS BViz wave- 
bands. We measure stellar masses for each pixel of each 
galaxy image from our BViz images by fitting to stellar 
population models. We use these values to construct stel- 
lar mass maps for our sample and compare morphologies 
in BViz and stellar mass. Our major findings and results 
include: 

I. We construct stellar mass maps and mass- to- light ra- 
tio maps for each galaxy and we present examples of each 
morphological type, in the zsso band, stellar mass maps, and 
[M/L)b in Figures [2]to 1101 We find that some compact el- 
liptical (cE) galaxies have blue cores, and more complicated 
internal structures than in either zsso or stellar mass would 
reveal on their own. Many of the peculiar ellipticals (pE) 
also have blue cores and we assert that these objects may 
have merged in their recent histories. The early and late- 
type spirals display more varied patterns. We find that for 
some of these galaxies, structures seen in light (e.g. spiral 
arms) are smoothed out in stellar mass. However, this effect 
does not hold true for all galaxies, nor for all features seen in 
light. Stellar mass and (M/L) maps of the 'peculiar' galax- 
ies are complicated and vary across the sample. Although 
it is more difficult to make out features in stellar mass for 
galaxies that have very blue colours, structures not seen in 
light are revealed. 

II. We compare the half-light radius, i?e, in zsso with 
half-mass radii for our sample as a function of morphologi- 
cal type. We find no systematic tendency for any particular 
morphological type to have larger Re in stellar mass than 
in the zsso band, and thus conclude that there is no bias 
introduced by measuring galaxy sizes in zsso- 

HI. We find a clear tendency for many galaxies to be 
more asymmetric in stellar mass than in zsso- We also find 
that morphology correlates with asymmetries, with early- 
types having low A(M*), and late- types higher values of 
A(M*). We find a relation between colour and A{M^) such 
that bluer galaxies have higher stellar mass asymmetry dif- 
ferences between optical light and stellar mass maps. The 
late- type spiral galaxies in the sample have higher A(M*) 
than would be expected from their asymmetries in ^sso • We 
discuss possible causes of this effect, including regions of en- 
hanced star formation also possessing higher stellar masses, 
and the evolution of spirals. We note that these highly 
asymmetric spirals resemble the clumpy disks of Elmegreen 
et al. (2007), and Conselice et al. (2004) and are experiencing 
either minor merging activity, or bulge formation through 
accretion of disk material. 

IV. We find that the Gini index in stellar mass is higher 
than in the z-band (G(M*) > G(^85o)) for all galaxies ex- 
cept the early- types, indicating that most of the stellar mass 
in later morphological types is contained within fewer pix- 
els. Late- type and peculiar morphologies show a trend for 
M2o(M*) > M2o(^85o), suggesting that the brightest 20 per- 
cent of the brightest pixels is not necessarily where the great- 
est 20 percent of the stellar mass is located. 

V. We investigate the relations between several combi- 
nations of the CAS parameters in BViz and stellar mass, 
and compare our results to previous morphological stud- 
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ies of this type (e.g. Conselice, Rajgor & Myers 2008; Lotz 
et al. 2004; 2008). We find tiiat 2:350 is the most appropriate 
photometric band to utihse for a z < 1 sample of galaxies 
with all morphologies, although stellar mass maps are bet- 
ter at distinguishing active galaxies from passive ones and 
is more a physical measure of structure. 

We furthermore compare our sample classifications in 
G-M20 to those of Lotz et al. (2008) and find that the Lotz 
et al. criteria do not best describe our sample. We revise 
the Lotz et al. (2008) criteria to best fit our sample in the 
zgso-band and within the stellar mass maps, and find that 
early-types, late-types and mergers can be separated. How- 
ever, the edge-on disk galaxies remain problematic and can- 
not be distinguished from mergers in G'(z85o)-M2o(^85o)- We 
find that G(M*)-M2o(M*) can be used to broadly separate 
early from late-type galaxies, but the criteria cannot distin- 
guish between late- type/edge-on disks from peculiar/merger 
systems. 

VL We find a relationship between and C (see ap- 
pendix A) for early-type galaxies in both 2^850 and stellar 
mass. In each case we find that galaxies with higher concen- 
trations have larger radii, and this relation is steeper in 2:550 
than stellar mass. We also investigate asymmetry versus M20 
and find that these parameters display a similar relation to 
A — C m. Z775 and 2:350- In stellar mass, A[M^)-y\.2o{M^) 
shows a tighter relation than A(M*) — C(M*), with a clear 
separation between early and late-type systems. Thus, we 
conclude that, for all parameters, late- type spiral galaxies 
overlap with with Pec/pM/M systems which may be ulti- 
mately taken from the same subset of galaxies. Structural 
studies in stellar mass do however track both minor and ma- 
jor merging events, and can be used to find galaxies in active 
galaxy evolution modes. 

We thank the GOODS team for making their data pub- 
lic, and STFC for a studentship towards supporting this 
work and the Leverhulme Trust for support. 
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APPENDIX A: APPENDICES 

APPENDIX B: STRUCTURE SCALING 
RELATIONS 

In this appendix we describe several additional scaling re- 
lationships between structural and morphological parame- 
ters, not discussed in the body of the paper, for our galaxy 
sample. Many of the relations we see here confirm the con- 
clusions we reached earlier. Furthermore, we list these for a 
complete overview of how structure correlates in the Hubble 
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Figure Bl. Asymmetry vs. dumpiness for the Vqqq, 2775 and zgbo bands. The sohd horizontal hne marks A = 0.35. The sohd 

diagonal line shows the fit for each band to the elliptical and spiral (i.e., 'normal') galaxies in the sample and the dashed line marks the 
position oi A = S. The black square point in the top right hand section of the plot illustrates the average error for the whole sample. 



Space Telescope ACS BViz bands as well as in stellar mass 
maps. 



Bl Asymmetry vs. dumpiness 

We postulate that the many small scale asymmetries, mea- 
sured by the dumpiness parameter, must approximate to 
the asymmetry of the galaxy as a whole. Figures IBll and 
IB2I show the relation between A and S for B435, Veoe, ^775, 
Z850 and stellar mass, respectively. The horizontal solid line 
in each plot marks the A — 0.35 position and the dashed 
line in each plot illustrates where A = S. The solid diagonal 
lines in each plot are the fits to the E, eS and IS galaxies 
in each band, which were forced to coincide at A = S = 0. 
These fits are shown below in Equations IB1[ IB5[ for each 
band respectively. 



^(^606) = 1.30(±0.08) • SiVeoe) 

A{i775) = 1.20(=t0.08) -5(^775) 

A{z85o) = 1.24(d=0.08) • S{z85o) 
A(M*) = 1.16(±0.06) • S{M.) 



(B2) 
(B3) 
(B4) 
(B5) 



A(B435) = 0.55(±0.07) • S{B435) 



(Bl) 



It is evident that the relations are mostly steeper than 
A — S, although the fit in stellar mass comes closest to this 
equality. The A{M^)-S{M^) relation thus is useful for iden- 
tifying merger candidates. The criterion of A{M^) > 0.35 is 
too low for the stellar mass maps, instead we suggest that 
galaxies which lie in the region above ^(M*) > 0.6 (horizon- 
tal dashed line. Fig. [B2} and underneath the A(M*)-5'(M*) 
fit, be classified as mergers, and this could be adapted for 
use to calculate merger fractions across a range of redshifts. 

The trend noted for the parameters to be less reliable in 
the bluer bands also holds here (Fig IB1|) . We note that the 
criteria A{z85o) > 0.35 and ^(2:850) > 1-24 • S{z85o) are ap- 
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Figure B3. The relation between C-M20 for the B435, Veoe, "^775 and zgbo bands. The sohd hnes show our fits to the E and eS galaxies 
in each band (Equations IBGl to IB9[) . The points are the same as in Figure Al. 



propriate for identifying merger candidates, and calculating 
merger fractions. The fits to A-S for ^775, zgbo and stellar 
mass are very similar, despite the much larger scatter in 
A(M*). 



B2 M20 vs. Concentration 

In this section we investigate the relationship between M20 
and C in B435, Veoe, ^775, zgbo and stellar mass (Figures IB3I 
and IB4p . The solid lines in each plot show the fits to the 
E and eS galaxies in the sample. The forms of these fits are 
shown below in Equations IB6l to IBIOI 

C{B435) = -0.92(±0.12) • M2o(B435) + 1.87(±0.16) (B6) 

aVeoe) = -1.32(±0.14) • M2o(V606) + 1.41(±0.20) (B7) 

C{i775) = -1.71(±0.13) • M2o(i775) + 0.44(±0.22) (B8) 

C7(z85o) = -1.48(=b0.13) • M2o(^85o) + 0.79(±0.22) (B9) 



C(M*) = -1.42(±0.08) • M2o(M*) + 0.81(±0.13) (BIO) 

M20 is anti-correlated with C, such that highly concen- 
trated galaxies have lower M20 values. The general trend 
for concentration to increase as M20 decreases holds for all 
bands, although the scatter increases at higher values of M20 
and lower C. The relation has a larger dispersion in the B435 
and Veoe bands, but this is more ordered than the random 
scatter seen in previous relations (Fig. IB3[ top panels). The 
relation is distinct in both 2775 and zsso , with the early- type 
galaxies occupying the high C, low M20 region, the eS's hav- 
ing higher M20 values for their concentration and late-types 
generally having low C and high M2o- 

This trend is also seen in the stellar mass maps (Fig. 
IB4|) , but with greater separation between the morphological 
types than in ^775 and zsso- The early- types also have high 
concentration and low M20, although the spirals are more 
randomly distributed, and do not occupy any one region 
of the plot. The late/Pec/pM/M galaxies have higher M20 
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S(Mass) 



Figure B2. Asymmetry vs. dumpiness for the stellar mass maps. 
The convention is the same as in Fig. IBll The horizontal dashed 
line marks the position of A = 0.6, which we derive as a criterion 
for merger candidates in stellar mass. 




M2Q(Mass) 

Figure B4. The relation between C-M20 for the stellar mass 
maps. The solid line shows our fit to the E and eS galaxies in 
stellar mass (Eq. IBlOf) . The solid line is this relation plus 2a. 
The points are the same as in Figure Al. 



values for their concentrations. The dashed line in Fig. IB4I 
shows the fit plus 2a and we postulate that galaxies falling 
above this criterion can be picked out as merger candidates. 



B3 Size vs. Concentration 

In this section we investigate the relation between galaxy 
size, as measured by the half-light radius. Re, and concen- 
tration, in ^850 and stellar mass. The top panel in Fig. IB5I 
shows the relation in zgbo and the bottom panel in stellar 
mass. We find that the early- types follow a linear relation- 
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Figure B5. The relation between Re and concentration for zgbo- 
band images (top) and stellar mass maps (bottom). The solid 
lines show the fits to the E galaxies in each band (Equations IBllI 
and IB12|) . The points are the same as in Figure Al. 



ship, and we have fit this relation for the E galaxies in zsso 

(Bll) 
(B12) 



(Eq. IBTTT) and stellar mass (Eq. [BT2]) . 
Re{z85o) = 0.78(=b0.15) • C(z85o) - 0.91(±0.52) 
Re(M,) = 0.48(±0.21) • C(M*) + 0.23(±0.68) 



The relation in zsso is steeper than in stellar mass, and both 
trends show that as Re increases, C also increases. For nor- 
mal galaxies, C is a tracer of stellar mass, such that lower 
C corresponds to less massive galaxies (Conselice 2003). In 
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Figure B7. Asymmetry-M2o for the stellar mass maps. The 
solid diagonal line marks our division between early and late- 
types/mergers. 

both ^850 and stellar mass we see that more concentrated 
galaxies have larger radii, and this effect is more pronounced 
in ^850 than in stellar mass. No such relation is found for the 
late- types and Pec/pM/M systems. 

Conselice & Arnold (2009) examine galaxies at z — 4 — 6 
in the Hubble Ultra Deep Field (HUDF) finding that early- 
types are smaller for their stellar mass at higher redshifts. 
They conclude that those galaxies that follow the C — Re 
relation are in a relaxed state. 

B4 Asymmetry vs. M20 

Here we investigate the relation between asymmetry and 
M20 in the B435, Veoe, ^775, 2:850 bands (Fig. IB6|) and stellar 
mass (Fig. IB7p . As M20 is analogous to C, we would expect 
to see a similar relation between A — M20 and A — C.ln Fig. 
IB6I we see that in the BViz bands we indeed find a similar 
relation in A — M20 to A — C, and could be used alongside 
A — C for classification purposes. 

In stellar mass, however, the relation between A — M20 
appears tighter than A{z85o) — C{z85o), such that galaxies 
with higher asymmetries have lower values of M20. There is 
a clear separation between early and late-type galaxies in 
stellar mass, and we have included an approximate divide in 
this plot (solid line in Fig. IB7p , the form of which is shown 
below in Equation IB13I 

A(M*) = -0.5 • M2o(M*) - 0.42 (B13) 

Although this divide is very successful in separating early- 
types from Pec/pM/M systems, the spirals are spread over 
these two regions, again implying that measuring structural 
parameters in stellar mass does not distinguish between mor- 
phologies as well, although it traces physical processes bet- 
ter. 

However, it is possible that viewing the stellar mass 
distribution of galaxies reveals different physical processes 
to those we expect to see in light. We have already discussed 



the possibility that some highly asymmetric late-type spirals 
in our sample could be in the process of forming bulges. It is 
also possible that studies in stellar mass could provide clues 
as to the formation mechanism involved. The asymmetric 
late- type spirals, for example object 19535 (Figs. [6l [TT]), 
could be undergoing outside-in formation, perhaps through 
the accretion of smaller satellites. We could, therefore, be 
picking out systems experiencing minor merging events. 



APPENDIX C: CHARACTERISATION OF 
DIFFERENT GALAXY TYPES 

At the start of this paper we describe how we classified each 
of our galaxies into one of nine type (§3.2). In this section, 
we described, based on our stellar mass map analysis, the 
properties of each of these galaxy types in terms of their 
optical and stellar mass structures, and what these imply 
about the evolution of different classes of galaxies. 

While galaxy classification in terms of overall visual 
morphology is limited (Conselice 2006b), it is still useful 
for connecting galaxies at different redshifts, and for under- 
standing how major nearby galaxy types such as disks and 
ellipticals are formed. Below we give details on the proper- 
ties of each of the nine galaxy type that we have examined 
within GOODS. 

CI Compact Es 

Compact ellipticals are typically found in the late-type area 
of the C-A optical light diagram, but are within the early- 
type region of the C-A stellar mass diagram. The dumpiness 
values for these systems are also zero in all wavelengths and 
stellar mass due to the fact that these systems are so small 
this parameter cannot be reliably measured. The concentra- 
tions of these systems are always quite low, with C = 2.5 — 3, 
at the resolution limit where we can measure these values, 
suggesting that they are not similar to giant ellipticals (e.g., 
Conselice 2003). 

Some interesting features appear when we compare the 
quantitative properties of cEs as measured in optical light 
to those in stellar mass maps. One of these is that for all 
cEs, the size measured in the stellar mass maps are always 
larger than for those measured in the reddest band, Z850- 
The asymmetry measured in stellar mass maps is also quite 
low, and lower than in the Z850 , perhaps due to the fact that 
the asymmetric light is brought to a lower level within the 
stellar mass images compared to the optical light. The M20 
parameter is also quite different between the stellar mass 
maps and the zgbo , such that M20 is higher (less negative) 
in the ^850 • This implies that the light distribution is more 
concentrated in the stellar mass maps, likely due to the blue 
coloured core. 

C2 Ellipticals 

First, we find that just based on colours, there is a wide 
diversity at a given redshift (see Figs. 12 and 13) - an in- 
dication that at higher redshifts there is a range of ellipti- 
cal galaxy star formation histories, something that has been 
recognised for some time (e.g., Stanford et al. 2004; Con- 
selice et al. 2007). Simply based on the C-A diagram with 
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Figure B6. Asymmetry-M2o for the B435, Veoe, "^TTS and 2:850 bands. The horizontal sohd hne marks A = 0.35 in each band. The points 
are the same as in Fig. Al. 



optical light (Fig. 20) it appears that most of the early- 
types are found in the region of C-A space classifiable as 
mid- types. This is commonly seen for high redshift ellipti- 
cals which appear more asymmetric than their z = coun- 
terparts (e.g., Conselice et al. 2007). Ellipticals do have gen- 
erally lower Gini and M20 indices, however there is still some 
spread in their values. 

The stellar mass maps reveal that ellipticals appear to 
become more 'diffuse' within stellar mass maps. However 
the sizes are either very similar in stellar mass maps or even 
slightly smaller, than in the corresponding 2;850-band image. 
The Gini and the M20 parameters do not change signifi- 
cantly between the stellar mass maps and the ^sso-band im- 
age. However, we find that some of the ellipticals are slightly 
more asymmetric in stellar mass maps than in the ^— band, 
although many of these systems tend to be bluer galaxies. 
The ellipticals also tend to become less concentrated in stel- 
lar mass compared with the zgbo band (Fig. I15p . 

Finally, the elliptical galaxies are one of the few types 
which do not change significantly in terms of their Gini and 
M20 values between stellar mass and the zsso, and this re- 



sults in these systems having a well defined location in the 
stellar mass Gini/M2o diagram (Fig. I23p . 



C3 Peculiar Ellipticals 

The peculiar ellipticals are a galaxy type which has been 
recognised to exist up to at least 2; = 1.4 before (e.g., Con- 
selice et al. 2007). These peculiar ellipticals appear elliptical 
in overall morphology, but have internal structure suggesting 
that they have recently undergone some type of assembly. 
These galaxies are not uncommon, and are more often seen 
than the compact elliptical type. 

These galaxies are more asymmetric and less concen- 
trated than the Es - with all having concentrations, C < 3. 
They are quite distinct from the ellipticals in the CA plane 
and are bluer by 0.5 magnitudes from the giant ellipticals. 
These blue features can be readily seen within the M/L maps 
for these systems (Fig. [3|). These galaxies often appear more 
elliptical in the stellar mass maps once their M/L variations 
have been taken into account. 

Another difference between the pEs and the Es is that 
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these systems have half-hght (or half-mass) radii which are 
nearly always bigger in the stellar mass maps than in the 
^850 band. This is a sign for star formation activity in the 
cores of these galaxies, since the bluer material in the core 
is less prominent in the stellar mass maps, giving them a 
larger effective radius. This can also be seen in the fact that 
the asymmetry in the stellar mass maps is generally higher 
than in the zgdo band for these systems, which is another 
indication that they are undergoing some form of star for- 
mation. 

The Gini values for these systems are higher in the stel- 
lar mass band than in zsso, which demonstrates that the 
stellar mass distribution is more concentrated in bright (but 
non-central) pixels than the light, as the light is less equally 
distributed than the stellar mass. The location of the pEs is 
similar to that of the Es in the Gini/M2o plane. 



types) do not change in size much between the stellar mass 
and the 2:850 image. However, the asymmetry in the stellar 
mass band is higher than for the ^350 band. This is likely due 
to the stellar mass within these galaxies being more concen- 
trated in fewer locations. The M20 also shows that the light 
is less concentrated in than the stellar mass, showing that 
the stellar mass is more distributed spatially than in light, 
suggestive of a merging systems. In fact, these peculiars have 
the highest asymmetries, M20, and Gini indices. 



C4 Spirals 

The spiral galaxies are amongst the most interesting for this 
study given that they often contain two major stellar pop- 
ulation types segregated spatially. Traditionally this is seen 
as an older stellar population making up the bulge or centre, 
and the spiral arms consisting of younger stellar populations. 

We find that the early- type spirals, those whose appar- 
ent bulge is brighter /larger than its apparent disk, are quite 
blue in colour. This difference is very obvious when examin- 
ing the M/L ratio maps for these galaxies (Fig. [5|), although 
the differences are not as obvious as with the late-type spi- 
rals. Nevertheless, these systems often appear to have a blue 
outer M/L ratio and a redder inner one, although there are 
examples where this is not the case. 

The sizes of these galaxies does not change significantly 
between the stellar mass maps and the zsso band images. 
There are however some differences between the stellar mass 
maps and the zgso band images within the various param- 
eters. The most significant of these is that the asymmetry 
parameter is higher in the stellar mass maps than in the 
^850 band imaging as discussed in detail in §4.3.3. The Gini 
parameter is also higher within the stellar mass maps than 
within the ^^sso by roughly 10 percent. This implies that the 
stellar mass is distributed in fewer pixels than in the light 
for spirals. The distribution of M20 values is similar between 
the light and stellar mass, suggesting that the spatial distri- 
bution is similar. 

One notable feature of the spiral galaxies in general is 
that they have a much larger range of asymmetry values 
within the stellar mass images. In fact, their values are of- 
ten as high as the peculiar galaxies, suggesting that within 
stellar mass maps these parameters cannot easily distinguish 
between mergers and spirals with mixed stellar populations. 



C5 Peculiars/Mergers 

The Peculiars, Pre-Mergers and Mergers all have similar pat- 
terns in the CAS space in both the stellar mass maps and 
within the 2:850 diagrams. These are the most asymmetric 
and bluest galaxies in our sample, and are likely in some 
phase of a merger, as has been described earlier in e.g., Con- 
selice et al. (2003, 2008). 

These peculiars (a term we use here for all three of these 



